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Introduction 

The purpose of the project is to refine and validate a probabilistic 
spatial ocrrputer model through the analyses of thematic mapper imagery. 

The model is designed to determine hew the interface between marshland and 
water changes as marshland is converted to wat$r in a disintegrating marsh. 
Coastal marshland in Louisiana is disintegratinr at the rate of approxima- 
tely 40 aq mi a year (Gagliano et al. 1981), and an evaluation of the 
potential inpact of this loss on the landings of estuarine-depe n dent 
fisheries is needed by fisheries managers. Understanding how marshland- 
water interface changes as coastal marshland is lost is essential to the 
prooess of evaluating fisheries effects, because several studies suggest 
that the production of estuarine-dependent fish and shellfish may be more 
closely related to the interface between marshland and water than to 
acreage of marshland (Faller 1979, Dow 1982, Zimmerman et al. 1984). The 
need to address this practical problem has provided an opportunity to apply 
sane scientifically interesting new techniques to the analyses of satellite 
imagery. Our progress with the development of these techniques is the sub- 
ject of our first report. 

The study group consists of two research teams. The first, located at 
the Southeast Fisheries Center in Miami (Miami Unit of the Beaufort, N. C., 
Laboratory) consists of Joan Browder and Alan Rosenthal. The second, 
located at Louisiana State University in Baton Rouge, La., consists of 
Nelson May, Robert Baixnann, and James Gosselink. Image analyses is being 
performed on the Fisheries Image Analyses System at the Slidell, La., 
Laboratory of the Mississippi Laboratories of the Southeast Fisheries 
Center. 
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Following are sumaries of work aooonpliahed to date and immediate 
future plana of the two research teams. 

Background 

The model simulates a disintegrating marsh, starting with a solidly 
vegetated r.tarsh and ending with total open water. Disintegration proceeds 
one pixel at a time, the specific pixel converting to water at each itera- 
tion determined by a probability function linked to a random raxrber genera- 
tor. Vteighting factors in the probability function allow us to wei^it the 
probability that a pixel will be converted at the next interation based on 
the nuiber of sides on vhich it is exposed to water or based on its posi- 
tion in the marsh. 

The probability function is: 

F-l+WS+GB 

where F ■ frequency of a given pixel on the selection list (relative 
probability) , W » weight given per side bordered by water, S • nrnber of 
sides bordered by water, G = weight given to pixels initially bordering the 
main water body, and B = a Boolean character equal to 1 or 0 that indicates 
whether a pixel borders the main water body. When W « 0 and G « 0, 
selection is randan except that boundary pixels have a slightly lower 
probability of disintegrating throughout the simulation because they are 
assured to have one side bordered by permanent land. F is calculated for 
each pixel at each iteration. 

Early modeling results (Browder et al. 1985) indicated that the 
marshland-water interface reaches its maxinun in an area vhen the spatial 
ocnpoeition of the marsh is about 50 percent land and 50 percent water but 
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that the magnitude of maximum interface changes inversely with the 
weighting factor of the probability function that converts segments of land 
to water. The strength of the weighting factor is reflected in the spatial 
distribution of land and water. With increasing weighting factor, the 
degree of clustering of water pixels is greater and water bodies become 
larger and less scattered. 

We plan to refine and validate our model by measuring and oaiparing 
the spatial distribution of land and water in simulated marshes and actual 
marshes in Louisiana, adjusting model coefficients so that spatial patterns 
in simulated marshes more nearly approximate those in the actual marshes, 
and placing actual marshes on disintegration curves (from zero to 100% con- 
version of land to water) produced from model simulations. The successful 
matching of model probability functions to patterns of land and water in 
actual marshes should not only quantify the relationship between magnitude 
of maxinun interface and patterns of land and water in the actual marshes 
but also confirm the model's prediction that maximum interface always 
occurs at about 50 percent land loss. 

A critical phase of the study is the development of quantitative para- 
meters that are descriptive of the spatial distribution of land and water. 
The development and testing of such parameters has been a major activity 
of the IWFS team for the past six months. Software development to measure 
both marshland-water interface and the parameters describing land and water 
parameters has been the major activity at LSU during its three months of 
involvement (since Septenfcer). 

fcWS Activities 

The behavior of the model was tested beyond the limits of the original 
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tests discussed in the preposed by extending the weighting factor for sides 
adjacent to water (W) up to 36 (originally tested it only through 5). 

We found that the stage in disintegration at Which maxinun interface occurs 
remained around 50 percent (Fig. 1, Table 1), but that maximum interface, 
expressed as percent maximum possible interface, decreased from about 51 
percent at W ■ 0 (the randan case) to 23 percent at W * 36 (Fig. 2, Table 
1 ). 

Algorithms were developed to measure spatial patterns. Four parame- 
ters vhich have been inp] eamented are: (1) frequency distribution of water 

pixels by water-cluster size; (2) frequency distribution of water clusters 
by size (in terms of water pixels); (3) frequency distribution of water 
pixels by number of sides bordered by other water pixels; and (4) percent 
pixels as water by distance frern main water body (in terms of pixel row). 
There are actually two versions of parameters one and two above: one 

including all water clusters and the other including only those water 
clusters connected to the main water body. 

In early tests with 20 repetitions of the same weighting factor, we 
found that parameter means appeared to reflect differences in the weighting 
factor. Means and confidence limits from these early tests are shewn in 
Figures 3 through 6. In these simulations, the marsh size vms 20 rows 
(pixels) by 50 oolums (pixels) . Pixels on each of the four edges of the 
marsh could have a maximum of three sides adjacent to water. Parameters 
were measured vrtien the marsh was 50 percent water. 

The frequency distribution of water pixels by water body size is shown 
for four different weiejiting factors in Figure 3. When W « 0, the randan 
case, more than half the water pixels are oontainsd in clusters of 1 to 25 
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water pixels. The nurber of water pixels contained in clusters of 1 to 25 
pixels decreases and the distribution moves to the ri^it as W increases. 

The frequency distribution of water pixels in terms of nurber of sides 
bordered by other water pixels is shown for weighting factors (W) from 0 to 
12 in Figure 5. The number of water pixels bordered on three or four sides 
by other water pixels increases with increasing weighting factor. 

The number of water pixels with distance from a man water body (in 
pixel lengths, or rows) is shown in Figure 5. With G set at 3, pixels in 
the first row (the row bordering the main water body) have a hi^ier proba- 
bility of disintegrating than pixels in the other rows. The result is that 
the front row has a higher nurber of water pixels at 50 percent disin- 
tegration than do the other rows. When W is increased first to 2 and then 
to 4, the second row also has a higher nurber of water pixels than the 
other rows. Otherwise, the nurber of water pixels in each row does not 
differ with distance. 

Work Scheduled for December 1985 to June 1986 

Further systematic tests of variation in spatial pattern will be 
conducted the spatial pattern parameters in future weeks. We are presently 
developing algorithms for measuring spatial autocorrelation. When 
carpi eted, the algorithms for measuring spatial patterns will be converted 
to FORTRAN 77 for inclusion in the ELAS software on tire FTPS system in 
Slidell. 

Our modeling work to date has been on a Hewlitt Packard 86-B rrdcroaoro- 
puter. We are rapidly reaching the memory limits of this ocrputer, even 
thou<b we have not yet expanded the size of the sixulated marsh, as we 
intend to do. In our recent work (beyond that covered ty this report), we 
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are interfacing the HP with a Burrou^is 6800 ndnioaiputer by tranefering 
files over a modem in order to increase execution sp eed. Our working mit 
eoa-i will receive an AT&T Unix-PC microocnputer, vAiich will be dedicated 
almost exclusively to this project and will give us greater capacity to 
expand both the size of the marsh and the number of parameters measured. 

To sunrrarize, activities in Miami during the next six months will con- 
sist of: 

* development of techniques for measuring spatial autocorrelation 

* measurement of spatial patterns on marshes simulated with alter- 
native weighting factors 

' conversion of spatial pattern parameters to FORTRAN 77 for inclu- 
sion in ELAS 

* conversion of the model program to C for execution on the AT&T Unix 
microcanputer 

* expansion of the size of the simulated marsh (in number of pixels). 

LSU Activities 

Work performed during the first three months of the project was con- 
centrated in two main areas: (1) activities associated with the start-155 

of the project, such as meetings with co-investigators and acquisition of 
magnetic tapes, ccnputer disk packs, topographic maps, and other supplies 
needed for the project, and (2) work with a contractor to begin some custom 
modifications of image-processing software required for the project. 

Image-prooessing for the project is being performed at the NASA 
Slidell Ccnputer Conplex, Slidell, Louisiana, vhere the Southeast Fisheries 
Center maintains the Fisheries Image Processing System (FIPS) . The system 
consists of a Sperry-Univac V77/600 oonputer and associated hardware and is 
equipped with a modified version of the Earth Resources Laboratory 


6 


Applcations Software (ELAS) (Graham «t al. 1984). ELAS is a FORTRAN-based 
image-processing package developed by NASA for processing digital data firm 
satellite remote sensors. Software modifications involved adapting 
existing capabilities in the latest version of ELAS to operate in the FIPS 
version of ELAS used in this study. 

Adaptation of the Water Body (WBOD) Module 

The water-body module (WBOD) is designed to categorize a user- 
specified class— for exartple, open water— into three classes s open bodies; 
small, closed bodies; and large, closed bodies. Information tabulated frcm 
repeated runs of the module will be used with other ELAS routines to deter- 
mine the size-frequency distribution of water bodies fron Landsat TM images 
of the Louisiana coastal marshes. Correct operation of the module was 
verified by running WBOD on a series of test files with known size- 
frequency distributions of water bodies. The FIPS version of WBOD became 
fully operational in mid-October. 

Adaptation of the Shoreline Length (SLIN) MDdule 

Adaptation of the shoreline length (SLIN) module to operate in FIPS 
has just been oonpleted. The SLIN module will be used to measure the 
length of the shoreline in land-water images derived frcm classifications 
of Landsat TW data. Algorithms that preceded the development of SLIN 
(Faller 1979) were not syntnetrical; that is, the reversed of land and water 
classes in a given image did not yield the sane length measurement (Dow and 
Pearson 1982). SUN uses a different algorithm to measure shoreline, 
avoiding the problem in the older software. 

The lengths of boundaries between two adjacent land-cover types, such 
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as land and water, are difficult to maasure accurately in raster images. 

SUN has a tende n cy to underestimate or overestimate the true shoreline 

length in a given area. The magnitude and direction of the bias is related 

to several factors (Dow and Pearson 1982) t 

sizes of the pixels ccnpoeing the map 

level of reticulation in the shoreline 

preprocessing techniques used to produce the 
land-water map frcm Landsat imagery. 

Thus, experimental techniques may have to be devised to correct for the 
biases in shoreline lengths derived from Landsat images of the Mississippi 
deltaic plain. 


Work Scheduled for Deoenber 1985 to June 1986 

delivery of the Landsat TM data tapes for the project 

conversion of the Landsat TM data to an ELAS format 

development of an experiment to quantify the magnitude and 
direction of the biases in shoreline-length measurements 
derived frcm the SUN module 

selection of sample sites in the Mississippi deltaic plain by 
marsh type (fresh, brackish, salt) and delta lobe, with consi- 
deration of cloud cover and quality of the TM imagery 

measurement of percentage area of open water, shoreline 
lengths, maxisun shoreline lengths, and determination of the 
distribution of water body sizes for each sanple rite 

anadyaia of U. S. Fish and Wildlife Service Habitat Maps 
(Wicker et ad. 1980) to determine historicad changes in land- 
water ratios and shoreline lengths between 1956 and 1978. 
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Table 1 


Means and 95 percent confidence limits of perce nt land loss at 
maxinun interface and maxinun interface expressed as percent 
maxinun possible interface for various vrei$iting factors *W) 
(CM3, m20). 


Percent Land Loss Maxinun Interface (%) 


w 

1CI 

Mean 

UC1 

LCI 

Mean 

UCI 

0 

49.40 

50.69 

51.98 

50.94 

51.39 

51.85 

1 

47.55 

49.16 

50.77 

41.41 

42.04 

42.66 

2 

46.20 

48.08 

49.95 

37.99 

38.71 

39.43 

3 

48.22 

50.15 

52.08 

35.92 

36.60 

37.29 

4 

45.57 

47.73 

49.83 

33.80 

34.45 

35.08 

5 

49.23 

51.36 

53.49 

33.00 

33.58 

34.15 

6 

48.49 

50.41 

52.32 

31.39 

31.92 

32.45 

7 

47.61 

49.94 

52.26 

30.79 

31.54 

32.30 

8 

46.61 

49.51 

52.41 

30.03 

30.60 

31.17 

9 

49.81 

52.04 

54.27 

29.45 

30.15 

30.85 

10 

45.75 

48.20 

50.65 

29.24 

29.86 

30.48 

11 

46.50 

48.74 

50.38 

28.57 

29.29 

30.02 

12 

47.08 

49.46 

51.84 

28.16 

28.76 

29.36 

16 

46.89 

49.17 

51.45 

26.56 

27.31 

28.07 

20 

48.17 

50.40 

52.62 

25.32 

25.97 

26.62 

24 

46.55 

49.61 

52.67 

24.91 

25.63 

26.33 

28 

47.86 

50.67 

53.48 

23.78 

24.47 

25.15 

32 

47.43 

49.88 

52.33 

23.30 

23. fi 

24.56 

36 

46.55 

50.74 

54.92 

22.34 

22.91 

23.47 
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1. Percent land converted to water at maximum interface versus 
weighting factor (W). Dot indicates man and vertical line 
indicates confidence range (r^20) . 

2 . Maxinun interface as percent maximum possible interface versus 
weighting factor (W). Dot indicates mean and vertical line 
indicates confidence range ( n*20 ) . 

3. Frequency distribution of water pixels by water cluster site 
(in tene of nunber of pixels) at 50 percent disintegration. 
Top of bar indicates mean and vertical line indicates con- 
fidence range (n*20). 

4. Frequency distribution of water pixels by nunber of sides bor- 
dered by other water pixels at 50 percent disintegration. 

5. Number of water pixels by distance from water (in pixel 
lengths, or rots) at 50 percent disintegration. 
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WEIGHTING FACTOR 

Figure 2. Maxinun interface as percent maximum possible interface versus 
weighting factor (W). Dot indicates mean and vertical line 
indicates confidence range (n=20). 
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Figure 3. Frequency distribution of water pixels by water cluster size 
(in terms of number of pixels) at 50 percent disintegration. 
Tap of bar i n di c a te s mean and vertical line indicates con- 
fidence range (nri2 0). 
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Figure 4. Frequency distribution of water pixels by ranter of sides bor- 
dered by other water pixels at 50 percent diu integration. 
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